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Abstract 



ff^ ' We present an updated global analysis of two-flavor MSW solutions to the 

^^ ■ solar neutrino problem. We perform a fit to the full data set corresponding to 

\^ , the 825-day Super-Kamiokande data sample as well as to Chlorine, GALLEX 

'nI" I and SAGE experiments. In our analysis we use all measured total event rates 



as well as all Super-Kamiokande data on the zenith angle dependence, energy 



0^ ' spectrum and seasonal variation of the events. We compare the quality of the 

^£j . solutions of the solar neutrino anomaly in terms of conversions of f e into active 

i-^ ' or sterile neutrinos. For the case of conversions into active neutrinos we find 

Q-i 
I ' that, although the data on the total event rates favours the Small Mixing 

tj'' Angle (SMA) solution, once the full data set is included both SMA and Large 

^S^ . Mixing Angle (LMA) solutions give an equally good fit to the data. We find 

that the best-fit points for the combined analysis are Am'^ = 3.6 x 10"^ eV^ 

K^ ; and sin^ 29 = 0.79 with xLn = 35.4/30 d. o. f and Am^ = 5.1 x lO'^ eV^ and 

;_( ' sin 26 = 5.5 x 10~^ with Xmin ~ 37.4/30 d. o. f. In contrast with the earlier 

504-day study of Bahcall-Krastev-Smirnov our results indicate that the LMA 

solution is not only allowed, but slightly preferred. On the other hand we 

show that seasonal effects, although small, may still reach 11% in the lower 

part of the LMA region, without conflict with the negative hints of a day-night 

variation (6% is due to the eccentricity of the Earth's orbit). In particular 

the best-fit LMA solution predicts a seasonal effect of 8.5%. For conversions 

into sterile neutrinos only the SMA solution is possible with best-fit point 

Am^ = 5.0 X 10"*^ eV^ and sin^ 26* = 3. x 10"^ and xLn = 40.2/30 d. o. f. 

We also consider departures of the Standard Solar Model (SSM) of Bahcall 

and Pinsonneault 1998 (BP98) by allowing arbitrary ^B and hep fluxes. These 

modifications do not alter significantly the oscillation parameters. The best fit 

is obtained for ^B/^BssM = 0.61 and hep/hepssM = 12 for the SMA solution 

both for conversions into active or sterile neutrinos and ^B/^BssM = 1-37 

and hep/hepssM = 38 for the LMA solution. 



I. INTRODUCTION 

It is already three decades since the first detection of solar neutrinos. It was realized from 
the very beginning that the observed rate at the Homestake experiment IQ was far lower 
than the theoretical expectation based on the standard solar model with the implicit 
assumption that neutrinos created in the solar interior reach the earth unchanged, i.e. they 
are massless and have only standard properties and interactions. In the first two decades of 
solar neutrino research, the problem consisted only of the discrepancy between theoretical 
expectations based upon solar model calculations and the observations of the capture rate 
in the chlorine solar neutrino experiment. This discrepancy led to a change in the original 
goal of using solar neutrinos to probe the properties of the solar interior towards the study 
of the properties of the neutrino itself. 

From the experimental point of view much progress has been done in recent years. 
We now have available the results of five experiments, the original Chlorine experiment 
at Homestake 0, the radio chemical Gallium experiments on pp neutrinos, GALLEX |^ 
and SAGE |P, and the water Cherenkov detectors Kamiokande [§ and Super-Kamiokande 
I^J^. The latter has been able not only to confirm the original detection of solar neutrinos at 
lower rates than predicted by standard solar models, but also to demonstrate directly that 
the neutrinos come from the sun by showing that recoil electrons are scattered in the direc- 
tion along the sun-earth axis. We now have good information on the time dependence of the 
event rates during the day and night, as well as a measurement of the recoil electron energy 
spectrum. After 825 days of operation, Super-Kamiokande has also presented preliminary 
results on the seasonal variation of the neutrino event rates, an issue which will become 
important in discriminating the MSW scenario from the possibility of neutrino oscillations 
in vacuum HJTO|. 

On the other hand there has been improvement on solar modelling and nuclear cross 
sections. For example, helioseismological observations have now established that diffusion is 
occurring and by now most solar models incorporate the effects of helium and heavy element 
diffusion [ll,!^- The quality of the experiments themselves and the robustness of the theory 



make us confident that in order to describe the data one must depart from the Standard 
Model (SM) of particle physics interactions, by endowing neutrinos with new properties. In 
theories beyond the Standard Model of particle physics neutrinos may naturally have exotic 
properties such as non-orthonormality [|13|, flavour-changing interactions [1^, transition 



magnetic moments ||T5| and neutrino decays |jT6|, the most generic is the existence of mass. 



While many of these may play a role in neutrino propagation and therefore in the explanation 
of the data |T^ it is undeniable that the most generic and popular explanation of the solar 
neutrino anomaly is in terms of neutrino masses and mixing leading to neutrino oscillations 
either in vacuum |l^ or via the matter-enhanced MSW mechanism [jl9 . 



In this paper we study the implications of the present data on solar neutrinos in the 
framework of the two-flavor MSW solutions to the solar neutrino problem. We perform a 
global fit to the full data set corresponding to 825 days of data of the Super-Kamiokande 
experiment as well as to Chlorine, GALLEX and SAGE. In our analysis we use as SSM 
the latest results from the most accurate calculation of neutrino fluxes by Bahcall and 



Pinsonneault ||20[ which incorporates the new normalization for the low energy cross section 



5*17 = 19^2 eV b indicated by the recent studies at the Institute of Nuclear Theory [^]. We 



have also considered the possibihty of departing from the SSM of BP98 by allowing a free 
normalization of the ^B flux and hep neutrino fluxes and we present the results we obtain 
when we treat these normalization as free parameters. 

We combined the measured total event rates at Chlorine, Gallium and Super- 
Kamiokande experiments with the Super-Kamiokande data on the zenith angle dependence, 
energy spectrum and seasonal variation of the events. The goal of such analysis is not only 
to compare the quality of the solutions to the solar neutrino anomaly in terms of flavour 
oscillations of z/g into active or sterile neutrinos but also to study the weight of different ob- 
servables on the determination of the underlying neutrino physics parameters as emphasized 
in Ref . |2^ . The outline of the paper is as follows. In Sec. |I| we present the basic elements 



that enter into our calculation of the observables and the definitions used in the statistical 
combination of the data. Section ^ contains our results of the allowed (or excluded) regions 
of oscillation parameters from the analyses of the different observables. The results on the 
allowed regions from the combined analysis of the total event rates is contained in Sec. p.llA. 



The constraints arising from the Super-Kamiokande searches for day-night variation of the 
event rates are discussed in Sec. [111B[ In Sec. p.11 C| we discuss the information which can 



be extracted from the distortion of the recoil electron energy spectrum measured by Super- 
Kamiokande. The restrictions arising from the preliminary Super-Kamiokande data on the 
seasonal variation of the event rates are studied in Sec. [111D| . In Sec. [111Il| we present our 
results from the global fit to the full data set and we determine the allowed range of oscil- 



lation parameters which are consistent with all the data. Finally in Sec. |^ we discuss the 
possible implications of our results for future investigations. 

Our results show that for oscillation into active neutrinos although the data on the total 
event rates favours the SMA angle solution, once the full data set is included both SMA and 
LMA give an equally good fit to the data. We find that the best-fit points for the combined 
analysis are Am^ = 3.6 x 10"^ eV^ and sin^ 26 = 0.79 with xLn = 35.4/30 d. o. f. and 
Am^ = 5.1 X lO"*^ eV^ and sin^ 29 = 5.5 x lO'^ with xLn = 37.4/30 d. o. f. We note that in 
contrast with the earlier 504-day study of Ref. ||22| our results indicate that the LMA solution 



is not only allowed, but actually slightly preferred. The existence of hints that the LMA 



MSW solution could be correct was also discussed in Ref. [^. On the other hand we find 
good quantitative agreement with the recent results of [S] as well as qualitative agreement 



with the old results of |]2^ based on a smaller sample. We show that seasonal effects may 
be no-negligible in the lower part of the LMA region, without conflict with the negative 
hints of a day-night variation. In particular the best-fit LMA solution predicts a seasonal 
effect of 8.5%, 6% of which is due to the eccentricity of the Earth's orbit. For conversions 
into sterile neutrinos we find that only the SMA solution is possible with best-fit point 
Am^ = 5.0 X 10"^ eV^ and sin^ 29 = 3.0 x 10-^ and xLn = 40.2/30 d. o. f. We also consider 
departures of the Standard Solar Model (SSM) of Bahcall and Pinsonneault 1998 (BP98) 
by allowing arbitrary ^B and hep fluxes. These modifications do not affect significantly 
the oscillation parameters. We find that the best fit is obtained for ^B/^Bssm = 0.61 and 
hep/hepssM = 12 for the SMA solution both for conversions into active or sterile neutrinos 
and ^B/^BssM = 1-37 and hep/hepssM = 38 for the LMA solution. 



II. DATA AND TECHNIQUES 

In order to study the possible values of neutrino masses and mixing for the MSW solution 
to the solar neutrino problem, we have used data on the total event rates measured at the 
Chlorine experiment in Homestake 0] , at the two Gallium experiments GALLEX and SAGE 
I^J^ and at the water Cerenkov detector Super-Kamiokande. Apart from total event rates 
we have in this case the zenith angle distribution of the events, the electron recoil energy 
spectrum and the seasonal distribution of events, all measured with their recent 825-day 
data sample p|]. 

We first describe our calculation of the different observables. For simplicity we consider 
the two-neutrino mixing case 

Ve = cos 6 ui + sin 61/2, Vx = — sin z/i + cos 9 z/2 , (1) 

where x can label either an active, x = /x, r, or sterile neutrino, x = s. In order to account 
for Earth regeneration effects, we have determined the solar neutrino survival probability 
Pee in the usual way, assuming that the neutrino state arriving at the Earth is an incoherent 
mixture of the vi and z/2 mass eigenstates. 

r> -pSun-pEarth , pSun pEarth fr\\ 

-' ee — -Tel -* le "r -' e2 -' 2e l^J 

where Pj^"" is the probability that a solar neutrino, that is created as i>e, leaves the Sun as 
a mass eigenstate Pi, and p^j^"^^^ is the probability that a neutrino which enters the Earth 
as Vi arrives at the detector as Ve- Similar definitions apply to P^^"" and Pg^"^*'^. 
The quantity P^^"" is given, after discarding the fast oscillating terms, as 

Pfr = 1 - i^et" = \ + ^\' PLz)cOs[2e,^{t,)] (3) 



where Plz denotes the improved Landau- Zener probability [^ and 6{to)m is the mixing 
angle in matter at the neutrino production point. In our calculations of the expected event 
rates we have averaged this probability with respect to the production point. The electron 
and neutron number density in the sun and the production point distribution were taken 



from Ref. 26 



In order to obtain the conversion probabilities in the Earth, P^f"^*'', we integrate the 
evolution equation in matter assuming a step function profile of the Earth matter density 
(the Earth as consisting of mantle and core of constant densities equal to the corresponding 
average densities, pm — 4.5 g/cm? and pc — 11.5 g/cm?). To convert from the mass density 
to electron and neutron number density we use the charge to nucleon ratio Z/A = 0.497 



for the mantle and Z/A = 0.467 for the core. In the notation of Ref. [^], we obtain for 

pEarth -i pEarth 

-* 2e ~ -■- -* le 

pEarth^^^ = (Zstney + {WiCosO + WssinOf (4) 

where 6 is the mixing angle in vacuum and the Earth matter effect is included in the formulas 
for Z, Wi and W3, which can be found in Ref. ||2^. p^j^^^^ depends on the amount of Earth 



matter travelled by the neutrino on its way to the detector, or, in other words, on its arrival 
direction which is usually parametrized in terms of the nadir angle, $, of the sun at the 
detector site. Due to this effect the survival probability is, in general, time dependent. This 
Earth regeneration effect is important in the study of the zenith angle distribution of events 
as well as in their seasonal variation ll9|,E8|. 



A. Rates 



Here we update previous analyses of solar neutrino data p2|,^ by including the recent 
825-day Super-Kamiokande data sample. We perform a dedicated analysis of the seasonal 
variation data and we anticipate its future role in discriminating between different solutions 
of the solar neutrino anomaly. Working in the context of the BP98 standard solar model of 
Ref. pO] we also allow for a free normalization of the ^B flux and of the hep flux. In our 



statistical treatment of the data we follow closely the analysis of Ref. ||29| with the updated 
uncertanties given in Refs. [0,^ . 

In our study we use the measured rates shown in Table |. For the combined fit we adopt 
the x^ definition: 



Xr 



«j 



1,3 



th 



Rr^'iRf 



Tjexp 



(5) 



where Rf^ is the theoretical prediction of the event rate in detector i and i?^^^ is the measured 
rate. The error matrix aij contains not only the theoretical uncertainties but also the 
experimental errors, both systematic and statistical. 

The general expression of the expected event rate in the presence of oscillations in ex- 
periment i is given by Rf^ : 



R: 



th 



E 0fe dE,\k{E,) X K,(E.)(Pee(^.,t)) 
A:=l,8 



(6) 



where Ey is the neutrino energy, (pk is the total neutrino flux and A^ is the neutrino energy 



spectrum (normalized to 1) from the solar nuclear reaction k [^| with the normalization 
given in Ref. ||20|. Here ae^i (o"x,j) is the z/g {vx-, x = fi, t) interaction cross section in the 



Standard Model |^ with the target corresponding to experiment i, and {Pee{Ey,t)) is the 
time-averaged z/g survival probability. 

For the Chlorine and Gallium experiments we use improved cross sections cTq, j(£') {a = 
e, x) from Ref. ^6 



For the Super-Kamiokande experiment we calculate the expected signal 



with the corrected cross section given in Sec. [II C 



The expected signal in the absence of oscillations, Rf^^^, can be obtained from Eq.(|^) by 
substituting Pee = 1. In Table | we also give the expected rates at the different experiments 
which we obtain using the fluxes of Ref. 



B. Day-Night Variation 



As already mentioned, in the MSW picture the expected event rates can be different 
when the neutrinos travel through the Earth due to the u^ regeneration effect |28|. As a 
result in certain regions of the oscillation parameters the expected event rates depend on 
the zenith angle of the sun as observed from the experiment site, since this determines the 
amount of Earth matter crossed by the neutrino on its way to the detector. 

The Super-Kamiokande Collaboration has studied the dependence of the event rates 
with the period of time along the day and the night. They present their results in the form 
of a zenith angle distribution of events. 



In our analysis we have used the experimental results from the Super-Kamiokande Col- 
laboration on the zenith angle distribution of events taken on 5 night periods and the day 
averaged value, shown in Table |IV| which we graphically reduced from Ref. ^ 

We define x^ fo^' the zenith angle data as: 

Xl = E ^ (7) 

where we have neglected the possible correlation between the errors of the different angular 
bins which could arise from systematical uncertainties. The factor a^ is included in order 
to avoid over-counting the data on the total event rate which is already included in Xr- 
We compute the expected event rate in the period i in the presence of oscillations as, 

1 rT{cOS^fnax i) f 

Rf = -^ ' drY. <Pu dE,Xk{E,) X [(Te,.(^.)(Pee(^., r)) (8) 

+cr,,,(E,)(l-(Pee(^.,r)))] 

where r measures the yearly averaged length of the period i normalized to 1, so Atj = 
T{cos(^rnax,i) — 'T{cos^min,i) = -500, .086, .091, .113, .111, .099 for the day and five night 
periods respectively. 

Super-Kamiokande has also presented their results on the day-night variation in the form 
of a day-night asymmetry. 

Ad IN = J^^~ ./^/ = -0.065 ± 0.031 (stat.) ± 0.013(syst.) (9) 

' Day + Night ^ ' ^ ^ ' ^ ' 

Since the information included in the zenith angle dependence already contains the day- 
night asymmetry, we have not added the asymmetry as an independent observable in our 
fit. Notice also, that being a ratio of event rates, the asymmetry is not a Gaussian-distributed 
observable and therefore should not be included in a x^ analysis. 



C. Recoil Electron Spectrum 

The Super-Kamiokande Collaboration has also measured the recoil electron energy spec- 
trum. In their published analysis |^2[ after 504 days of operation they present their results 



for energies above 6.5 MeV using the Low Energy (LE) analysis in which the recoil energy 
spectrum is divided into 16 bins, 15 bins of 0.5 MeV energy width and the last bin containing 
all events with energy in the range 14 MeV to 20 MeV. Below 6.5 MeV the background of 
the LE analysis increases very fast as the energy decreases. Super-Kamiokande has designed 
a new Super Low Energy (SLE) analysis in order to reject this background more efficiently 
so as to be able to lower their threshold down to 5.5 MeV. In their 825-day data |p they 
have used the SLE method and they present results for two additional bins with energies 
between 5.5 MeV and 6.5 MeV. 

In our study we use the experimental results from the Super-Kamiokande Collaboration 
on the recoil electron spectrum on the 18 energy bins including the results from the LE 

6 



analysis for the 16 bins above 6.5 MeV and the results from the SLE analysis for the two 
low energy bins below 6.5 MeV, shown in Table |Tl|. 

Notice that in Table ^ we have symmetrized the errors to be included in our y^ analysis. 
We have explicitly checked that the exclusion region is very insensitive to this symmetriza- 
tion. We define x^ for the spectrum as 



' pBP98 « ' ^] v"«P pBP98 



where 



Again, we introduce a normalization factor a^p in order to avoid double-counting with the 
data on the total event rate which is already included in x|.. Notice that in our definition 
of x| we introduce the correlations amongst the different systematic errors in the form of a 
non-diagonal error matrix in analogy to our previous analysis of the total rates. These cor- 
relations take into account the systematic uncertainties related to the absolute energy scale 
and energy resolution, which were not yet available at the time the analysis of Ref. |^ was 
performed. Note that our procedure is different from that used by the Super-Kamiokande 
collaboration. However we will see in Sec. [111C| that both methods give very similar results 



for the exclusion regions. This provides a good test of the robustness of the results of the 
fits, in the sense that they do not depend on the details of the statistical analysis. 

The general expression of the expected rate in the presence of oscillations E!'^ in a bin, 
is given from Eq.(|^) but integrating within the corresponding electron recoil energy bin 
and taking into account that the finite energy resolution implies that the measured kinetic 
energy T of the scattered electron is distributed around the true kinetic energy T' according 
to a resolution function Res(T, T') of the form [p^ ]: 



Res(T, T') = — exp 

v27rs 



^ ^ (12) 



2s2 



where 



s = soVT'/MeV , (13) 

and So = 0.47 MeV for Super-Kamiokande 0,0] • On the other hand, the distribution of the 
true kinetic energy T' for an interacting neutrino of energy E^y is dictated by the differential 



cross section daa{E^, T')/dT', that we take from |^. The kinematic limits are: 



o<r<T'(£j, r'(E„) = ^^^. (14) 

For assigned values of Sq, T^in, and Tmax, the corrected cross section a^iE) (a = e, x) is 
given as 

/•Tniax fT'{E,,) da (E T') 

aaiE,) = dT dT' ResiT, T') ^^"7;; ^ . (15) 

JT,^i„ Jo dl 

In Fig. ^ we show our results for the recoil electron spectrum in the absence of oscillations 
and compare it with the expectations from the Super-Kamiokande Monte-Carlo. We see 
that the agreement is excellent. In this figure no normalization has been included. 



D. Seasonal Variation 

After 825 days of operation, Super-Kamiokande has also presented preliminary results 
on the seasonal variation of the neutrino event rates P] which seem to hint at a seasonal 
variation of the data, especially for recoil electron energies above 11.5 MeV (see Table |rv| ). 
As discussed in Ref. P,p!0| the expected MSW event rates do exhibit a seasonal effect. 
In the LMA region such dependence can be expected mainly due to the different night 
duration throughout the year at the experimental site and also due to the different averaged 
Earth densities crossed by the neutrino during the night periods, which lead to a seasonal- 
dependent z/g regeneration effect in the Earth. On the other hand, in the SMA region. Earth 
matter effects are due to resonant conversion of neutrinos in the Earth which is only possible 



when neutrinos travel both through the mantle and the core [^. Thus we find that in the 
SMA region the seasonal variation is associated with the fact that at Super-Kamiokande site 
only in October-March nights there are neutrinos arriving at sufficiently low zenith angle to 
satisfy the resonant condition. 

We define x^ for the seasonal variation data as, 

Xlea = E ^2 (16) 

j=l,8 '-'i 

where, as before, the normalization factor asea is introduced to avoid double-counting. 

Taking into account the relative position of the Super-Kamiokande setup in each period 
of the year, we calculate the distribution of the events as 

rU+At/2 ^^ 

Here At = 1.5 months and R^^{t) is obtained from Eq.(^ but using the time-dependent 
survival probabilities and integrating the recoil electron energy above 11.5 MeV. Notice that, 
unlike in the 708 days data sample, in order to compare our results with the recent data on 
the seasonal dependence of the event rates from the Super-Kamiokande Collaboration for 
the 825 data sample, we have not included the geometrical seasonal neutrino flux variation 
due to the variation of the Sun-Earth distance arising from the Earth's orbit eccentricity 
because the new Super-Kamiokande data is already corrected for this geometrical variation. 

III. FITS: RESULTS 

We now turn to the results of our fits with the observables described above. We have ob- 
tained the regions of allowed oscillation parameters Am^-sin^ 29 by obtaining the minimum 
X^ and imposing the condition x^ < Xmm+Ax^(2, CL) where Ax^(2,CL) = 4.61(9.21) for 
90% (99%) CL regions. 



A. Rates 

We first determine the allowed range of oscillation parameters using only the total event 
rates of the Chlorine, Gallium and Super-Kamiokande experiments. The average event rates 
for these experiments are summarized in Table |. We have not included in our analysis the 
Kamiokande data |^ as it is well in agreement with the results from the Super-Kamiokande 
experiment and the precision of this last one is much higher Q. For the Gallium experiments 
we have used the weighted average of the results from GALLEX |Q and SAGE ||^ detectors. 

Using the predicted fluxes from the BP98 model the x^ fo^' the total event rates is 
XssM — 62.4 for 3 d. o. f. This means that the SSM together with the SM of particle 
interactions can explain the observed data with a probability lower than 10~^^! 

In the case of active-active neutrino oscillations we find that the best-fit point is obtained 
for the SMA solution with 



i^ = 5.6 X 10-^eV^ 

,2 na « Q XX 1 n— 3 



sin^ 2^ = 6.3 X 10~% (18) 

2 f\ 07 

X,min ' 

what implies that the solution is acceptable with a 55% CL. 

There are two more local minima of x^- One is the LMA solution with the best fit for 

Am^ = 1.4 X IQ-^eV^ , 
sin^ 2^ = 0.67, (19) 

v^ = 2 92 

which is acceptable with at 91% CL. The other is the LOW solution with best-fit point 

Am^ = 1.3 X IQ-^eV^ , 
sin^ 29 = 0.94 , (20) 

y^ • =74 

which is only acceptable at 99% CL. 

In the case of active-sterile neutrino oscillations the best-fit point is obtained for the 
SMA solution with 



(21) 



acceptable at 90% CL. The LMA and LOW solutions are not acceptable for oscillation 
into sterile neutrinos. In those regions Xmm > 19.5 implying that they are excluded at the 
99.999% CL. Unlike active neutrinos which lead to events in the Super-Kamiokande detector 
by interacting via neutral current with the electrons, sterile neutrinos do not contribute to the 
Super-Kamiokande event rates. Therefore a larger survival probability for ^B neutrinos is 
needed to accommodate the measured rate. As a consequence a larger contribution from ^B 
neutrinos to the Chlorine and Gallium experiments is expected, so that the small measured 

9 



Am^ 


= 5.0x 


10' 


-^eV^ 


sin^ 29 


= 5.0x 


10" 


-3 


A-m«n 


= 2.6, 







rate in Chlorine can only be accommodated if no '^ Be neutrinos are present in the flux. 
This is only possible in the SMA solution region, since in the LMA and LOW regions the 
suppression of ^ Be neutrinos is not enough. 

In Fig. ^ we show the 90% and 99% CL allowed regions in the plane A?7i^-sin^ 26. The 
best-fit points in each region are marked. We find that as far as the analysis of the total 
rates is concerned, there is no substantial change in the best-fit points in the three regions as 
compared to the previous most recent analysis including the 504 days of Super-Kamiokande 
data ig. 



We have also considered the possibility of departing from the SSM of BP98 by allowing 
a free normalization of the ^B flux and we treat this normalization as a free parameter /? 
in our analysis. Figure |^ shows the allowed regions in the MSW parameter space when j3 is 
allowed to take arbitrary values. The best fit SMA solution is obtained for 

Am^ = 5.0 X 10-^eV^ , 
sin^ 29 = 5.0 x 10"^ , (22) 

P = 0.82 

y2 . = 0.05 . 

The best fit for the LMA solution occurs at: 

Am^ = L6 X 10"^eV^ , 
sin^ 26 = 0.63 , (23) 

p = L32 
Y^ • = 47 

and the LOW solution has its best-fit point at /3 = .98 and therefore coincides with the one 
obtained in Eq.(pO|). 

In the case of active-sterile neutrino oscillations the best-fit point is obtained for the 
SMA solution : 

Am^ = 5.0 X 10~^eV^ , 
sin2 2^ = 3.2 X 10-^ (24) 

p = 0.75 
Y^ • =2 16 

For all solutions we find that the main effect of considering a free normalization of the 
^B flux is upon the quality of the fits, as measured by the depth of the x^- Next comes the 
position of the best-fit points, mainly a reduction in the value of the mixing angle. The 
allowed regions are considerably enlarged as can be seen by comparing Figs. |^ and ^ 

B. Yearly Averaged Zenith Angle Dependence 

We now study the constraints on the oscillation parameters from the Super-Kamiokande 
Collaboration measurement of the zenith angle distribution of events. We use the data taken 
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on 5 night periods and the day averaged value shown in Table |I^ which we graphically 
reduced from Ref. ^. 

Considering only the zenith angle data (not including the information from the total 
rates), the best-fit point in the case of active-active neutrino oscillations is obtained for 
Am^ = 4.5 X 10"^ eV^ and sin^26' = 1.0 with xLn = 2.3 for 3 d. o. f. and in the case of 
active-sterile neutrino oscillations is obtained for Am"^ = 3.2 x 10~^ eV^ and siia^ 29 = .98 
with Xmin = 2.2. With these values we calculate the excluded region of parameters at the 
99% CL, shown in Fig ^. Notice that the zenith angle data favours the LMA solution of the 
solar neutrino problem. 

When combining the information from both total rates and zenith angle data, we obtain 
that in the case of active-active neutrino oscillations the best-fit point is still obtained for 
the SMA solution with 



Am^-- 


= 5.0x 


10- 


-^eV^ 


sin^ 29 -- 


= 6.3x 


10^ 


-3 




= 5.9, 







(25) 

for 6 d. o. f. which is acceptable with a 56% CL. However the LMA solution becomes 
relatively better with a local minimum at 

Am^ = 4.5 X lO'^eV^ , 
sin^ 2^ = 0.80, (26) 

Y^ =72 

valid at 70% CL. This difference arises from the fact that, although small, some effect is 
observed in the zenith angle dependence which points towards a larger event rate during the 
night than during the day, and that this difference is constant during the night as expected 



for the LMA solution ES]. In the SMA solution, however, the enhancement is expected to 



occur mainly in the fifth night [pS 



The LOW solution is almost un-modified and presents the best-fit point at 



Am^ = 1.0 X lO^^eV^ , 
sin^ 29 = 0.94 , (27) 

v^ • =12 7 

A.mm ' ' 

which is acceptable at 95% CL. 

In the case of active-sterile neutrino oscillations the best-fit point is obtained for the 
SMA solution with 



Am2 = 


= 5.0 


X 


10" 


-^eV^ , 


sin^ 29 -- 


= 5.0 


X 


10" 


~3 

5 


Amjn 


= 8.1 


) 







(28) 



valid at 77% CL. 

Figure ^ shows the regions excluded at 99% CL by the zenith angle data alone, together 
with the 90% and 99% CL allowed regions in the plane Am^-sin^ 29 from the combined 
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analysis of rates plus zenith angle data. The best-fit points in each region are indicated. 
The main difference with respect to Fig. ^ is observed in the LMA allowed region which is 
cut from below, as the expected day-night variation is too large for smaller neutrino mass 
differences. 



C. Recoil Electron Spectrum 

We now present the results of the study of the recoil electron spectrum data observed in 
Super-Kamiokande. Using the method described in Sec. |IIC| we obtain that the x^ for the 



undistorted energy spectrum (Standard Model case) is 20.1 for 17 d. o. f.. This corresponds 
to an agreement of the measured with the expected energy shape at the 27% CL. This 
value depends on the degree of correlation allowed between the different errors. In this 
way, reducing the correlation in the error matrix the agreement decreases down to 13%. 
This suggests that the correlations amongst the different systematic uncertainties related 
to the energy resolution will play an important role in the analysis of the energy spectrum. 
Indeed in our definition of Xs ^^ introduce the correlations amongst the different systematic 
errors in the form of a non-diagonal error matrix in analogy to our previous analysis of the 
total rates. These correlations take into account the systematic uncertainties related to 
the energy resolution. We believe that these might constitute the main difference between 



our treatment of the spectrum data and the earlier one presented in Ref. [^, when this 
experimental information was still unavailable. 

We find that the best fit to the spectrum in the MSW plane for the case of active- 
active neutrino oscillations is obtained for Am"^ = 6.3 x 10^^ eV^ and sin^ 29 = 0.08 with 
Xmin — 17.9. For the case of active-sterile neutrino oscillations we get Am^ = 6.3 x 10~^ 
eV^ and sin^ 29 = .08 with xLin = I'''- With these values we obtain the excluded region of 
parameters at the 99% CL shown in Fig. ^. 

We see that our results for the exclusion regions are quantitatively very similar to those of 
the Super-Kamiokande collaboration , even though our procedure is different from theirs. 
This provides a good test of the robustness of the results of the fit of the spectrum, in the 
sense that they do not depend on the details of the statistical analysis. In contrast we note 



that our results are different from those of Ref. |2^ . 

When combining the information from both total rates and the recoil energy spectrum 
data we obtain that in the case of active-active neutrino oscillations both SMA and LMA 
solutions lead to fits to the data if similar quality. In this way, the best-fit for the LMA 
solution is 



(29) 



while for the SMA solution we find 



Am^ 


= 1.4 X 


10 


-^eV^ , 


sin^ 29 


= 0.67 


J 




Amin 


= 22.5 


•) 




Am2 = 


= 5.6x 


10" 


6eV^ 


sin^ 29 = 


= 5.0x 


10" 


3 


A. mm 


= 23.4 , 







(30) 
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Am2 = 


= 5.0x 


10" 


-^eV^ , 


sin^ 29 -- 


= 3.0x 


10" 


-3 


A mm 


= 26.3 . 







for 18. d. o. f. which are acceptable at 83%. Finally for the LOW solution we find 

Am^ = 1.0 X lO^^eV^ , 
sin^ 29 = 0.94 , (31) 

v^ • = 26 7 

Arrwn ^' ' ' 

with agreement with data only at the 8.5% CL. 

In the case of active-sterile neutrino oscillations the best-fit point is obtained for the 
SMA solution with 



(32) 



acceptable at 90%. 

In Fig. H we plot the excluded region at 99% CL by the energy spectrum data together 
with the 90% and 99% CL allowed regions in the plane Am^-sin^ 29 from the combined 
analysis. The best-fit points in each region are marked. 

The main point here is that the oscillation hypothesis does not improve considerably the 
fit to the energy spectrum as compared to the no-oscillation hypothesis. In this connection 
it has been suggested [^] that better description can be obtained by allowing a larger flux 
of hep neutrinos as they contribute mainly to the end part of the spectrum. In order to 
account for this possibility we have also analysed the data allowing for a free normalization 
of the ^B and hep fluxes, treating them as a free parameters (3 and 7 correspondingly. When 
doing so we find that the no-oscillation hypothesis gives x^ = 17.4 for 15 d. o. f. for j3 = 0.45 
and 7 = 13.5. 

When combining the information from both total rates and the recoil energy spectrum 
data in the case of active-active neutrino oscillations we obtain that for the LMA: 

Am^ = 1.6 X lO'^eV^ , 
sin^ 29 = 0.63 , (33) 

/? = 1.3 7 = 33 
Y^ • = 17 5 

for 16 d. o. f. which is acceptable at 66% while for the SMA solution 

Am^ = 5.0 X lO^W^ , 
sin^ 29 = 2.5 x 10"^ , (34) 

/5 = .61 7 = 13 

y2 . = 19 5 

which is acceptable at 75% CL. Finally for the LOW solution we find 

Am^ = 1.0 X 10"^^ , 
sin2 2^ = 0.94, (35) 

/3 = .97 7 = 22 
v^ = 25 

Amm " ? 
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acceptable at 93%. 

In the case of active-sterile neutrino oscillations the best -fit point is obtained for the 
SMA solution : 

Am^ = 5.0 X lO^^eV^ , 
sin2 2e = 2.0 X 10"^ (36) 

(3 = 0.61 7 = 12 
y^ • = 22 . 

which is acceptable at 89% CL. 

In Fig. ^ we display the normalized expected energy spectra for SMA, LMA solutions 
for active-active oscillations and for no-oscillation with non-standard ^B and hep fluxes. 

D. Seasonal Dependence 

Recently Super-Kamiokande has also presented preliminary results which seem to hint 
for a seasonal variation of the event rates, especially for recoil electron energy above 11.5 
MeV. As explained in Sec. |11D| the expected MSW event rates do exhibit a seasonal effect 
due to the the seasonal- dependent v^ regeneration effect in the Earth. We explore here the 
constraints on the MSW oscillation parameters which can be extracted from the seasonal 
dependence data given in Table fV[ 

Considering only the seasonal variation data above 11.5 MeV (allowing a free normal- 
ization for the ^B flux) we find that the the SSM yields a value of x^ = §■ for 7 d. o. f. This 
shows that the data is still not precise enough to enable one to draw any definitive conclu- 
sion. However one can still obtain some preliminary information on the MSW parameters 
from the analysis of these data. In this way, when allowing oscillations into active flavours 
we obtain the best-fit point for Am^ = 3.2 x 10~^ eV^ and sin^26' = .1 with Xmin = 1-8 
for 5 d. o. f. and in the case of active-sterile neutrino oscillations best fit is obtained for 
Am^ = 1.3 X 10"*^ eV2 and sin2 2^ = 0.1 with xLn = 2-0. 

In Fig. 1^ we show the exclusion regions at 95% and 99% CL. The larger area represents 
the allowed region at 95% CL. Being an allowed region, it means that at 95% CL, some 
effect is observed. The darker area shows the small excluded region at 99% CL. 

Since the seasonal variation of the event rates in the MSW region is due to neutrino 
regeneration in the Earth one expects a correlation with the day-night variation which arises 
from the same origin. This correlation is observed as the 95% allowed regions in Fig. ^ have 
a large overlap with the 99% excluded region from the observed zenith angle dependence in 
Fig. ^. Notice however that the upper part of the 95% CL allowed region for the oscillation in 
active flavours (larger Am^ and larger mixing angles) is still in agreement with the observed 
zenith angle dependence. Thus, should the seasonal effect be confirmed in the higher energy 
part of the spectrum, it will favour the LMA solution to the solar neutrino problem. 

E. Combined Analysis 

We now present our results for the simultaneous fits to all the available data and ob- 
servables. In the combination we define the global x^ as the sum of the different x^ defined 
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above. In principle such analysis should be taken with a grain of salt as these pieces of in- 
formation are not fully independent; in fact, they are just different projections of the double 
differential spectrum of events as a function of time and energy. Thus in our combination 
we are neglecting possible correlations between the uncertainties in the energy and time 
dependence of the event rates. 

From the full data sample we obtain that the for the SSM of Ref. [^ 

xLniSSM) = 96. (37) 

for 32 d. o. f. So the probability of explaining the full data sample as a statistical fluctuation 
of the the SSM together with the SM of particle interactions is smaller than 10~^. 

In the MSW oscillation region we obtain that for oscillations into active flavours the best 
fits are obtained for the LMA solution with 

Am^ = 3.65 X lO^^eV^ , 
sin2 2^ = 0.79, (38) 

v^ = S5 4- 

Amin ^^-^ , 

for 30 d. o. f. which implies that the solution is acceptable at 77% CL., while in the SMA 
region the local best-fit point is 

Am^ = 5.1 X 10~^eV^ , 
sin^ 29 = 5.5 x 10"^ , (39) 

y^ • = 37 4 

valid at the 83% CL.. The global analysis still presents a minimum in the LOW region with 
best-fit point: 

Am^ = LO X IQ-W^ , 
sin^ 29 = 0.94 , (40) 

Y^ • =40 

which is acceptable at 90% CL. 

The results of the global analysis for the case of of active-sterile neutrino oscillations, 
show the fit is slightly worse in this case than in the active-active oscillation scenario. The 
best -fit point is obtained for the SMA solution with 

Am^ = 5.0 X 10"^^ , 
sin2 2^ = 3.0 X 10"^ (41) 

v^ ■ = 40 2 

acceptable at 90% CL. 

In Fig. H we show the 90% (lighter) and 99% (darker) CL allowed regions in the plane 
Am^-sin^ 29. Best-fit points in each regions are also indicated. 

By comparing Figs. ^ and ^ we see the effect of the inclusion of the full data from 
Super-Kamiokande on both time and energy dependence of the event rates. For oscillations 
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into active neutrinos, the larger modification is in the LMA solution region which has its 
lower part cut by the day-night variation data while the upper part is suppressed by the 
data on the recoiled energy spectrum. The best fit point is also shifted towards a larger mass 
difference by a more than a factor 2. As for the SMA solution region the position of the 
best fit point is shifted towards a slightly smaller mixing angle, while the size of the region 
at the 90% CL is reduced. At 99% CL the allowed SMA region is very little modified. For 
oscillations into sterile neutrinos the best fit point is also shifted towards a smaller mixing 
angle. 

We finally study the allowed parameter space with the normalization of the ^B and hep 
fluxes are left free. Figure |^ shows the allowed regions in the MSW parameter space when 
the normalizations are allowed to take arbitrary values. The best fit for the LMA solution 
occurs at 

Am^ = 3.6 X lO^^eV^ , 
sin^ 29 = 0.67 , (42) 

p = 1.37 7 = 38 

which is acceptable at 64% CL. The best fit SMA solution is obtained for 

Am^ = 5.0 X lO^^eV^ , 
sin^ 29 = 2.5 x 10"^ , (43) 

P = 0.61 7 = 12 
y' • = 34 . 

which is acceptable at 80% CL., The best fit for the LOW solution occurs at 

Am^ = 1.0 X 10~^eV^ , 
sin^ 29 = 0.94 , (44) 

/3 = 0.97 7 = 21 

v^ = "^H 1 
A.min oo.L , 

which is acceptable at 90% CL., 

In the case of active-sterile neutrino oscillations 

Am^ = 5.0 X 10"^eV^ , 
sin^ 29 = 2.0 x 10"^ , (45) 

p = 0.61 7 = 12 
v^ = 35 8 

A.min ^^.'-j . 

which is acceptable at 85% CL., 

For all solutions we find that the main effect of considering free normalization of the ^B 
and hep fluxes is upon the position of the best-fit points, mainly a reduction in the value of 
the mixing angle. The allowed regions are also enlarged as can be seen by comparing Figs. ^ 
and|. 
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IV. SUMMARY AND DISCUSSION 

We have presented an updated global analysis of two-flavor MSW solutions to the solar 
neutrino problem using the full data set corresponding to the 825-day Super-Kamiokande 
sample plus Chlorine, GALLEX and SAGE experiments. In addition to all measured total 
event rates we included all Super-Kamiokande data on the zenith angle dependence, energy 
spectrum and seasonal variation of the events. We have given a comparison of the quality 
of different solutions of the solar neutrino anomaly in terms of MSW conversions of z/g into 
active and sterile neutrinos. For the case of conversions into active neutrinos we have found 
that once the full data set is included both SMA and LMA solutions give an equally good fit 
to the data. We find that the best-fit points for the combined analysis are Am'^ = 3.6 x 10~^ 
eV2 and sin2 20 = 0.79 with xLn = 35.4/30 d. o. f and Am^ = 5.1 x lO'^ eV^ and 
sin^ 29 = 5.5 x 10^^ with Xmin = 37.4/30 d. o. f. In contrast with the earlier 504-day study 
of Bahcall-Krastev-Smirnov our results indicate that the LMA solution is slightly preferred. 
Although small, there is a hint for seasonality in the data and this should therefore take 
into account in future studies, as we have indicated here. Defining the seasonal variation (in 
percent) as 

jDth T?^^ 



R-^ + R 



th 

"min 



where Rmax(min) is ^^^ expected event rate during the winter (summer) period, we find 
that the seasonal effect may still reach 11% in the lower part of the allowed LMA region 
(6% comes from the variation expected from the geometric effect due to the eccentricity of 
the Earth's orbit), without conflict with the negative hints of a day-night variation. Such 
seasonal dependence is correlated with the day-night effect and this in turn can be used in 
order to discriminate the MSW from the vacuum oscillation solution to the solar neutrino 
anomaly. We have performed a numerical study of this correlation |T^ which generalizes 
the estimate presented by Smirnov at the 1999 edition of Moriond for a constant Earth 
density. Our results show that due to the Earth matter profile the seasonal variation can 
be substantially enhanced or suppressed as compared to the expected value obtained with 
an average Earth density. One must notice, however, that the exact values of masses and 
mixing for which a significant enhancement is possible may depend on the precise model 
of the Earth density profile. The numerical results presented above were obtained using 
the step function profile and thus numerical differences with the results obtained with, for 



instance the PREM model |23| can be expected for a given point in the MSW plane. 

With the results from the combined analysis one can also predict the expected event 
rates at future experiments. For example, we find that the average oscillation probabilities 
for ^Be neutrinos in the different allowed regions for MSW oscillations into active neutrinos 
is 

PrBe SMA = 0.01+°:J^4, (46) 

RBe LMA = 0.55+°:°t (47) 

These results imply that, at Borexino we expect a suppression on the number of events as 
compared with the predictions of the SSM of 
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^«MA/^BF98 ^ 0.22+^:^^4, (48) 

^LMA/^BP98 ^ o.eS^He- (49) 

For conversions into sterile neutrinos only the SMA solution is possible with best-fit point 
Am^ = 5.0 X 10-6 eV^ and sin^ 26 = 3.2 x lO^^ and xLn = 40.2/30 d. o. f which leads to 

-P^Bc sterile = 0.015_o!o02) (^O) 

and one expects a larger suppression of events at Borexino 

^sterile/^BP98 ^ 0.015+°:°^0. (51) 

As a way to improve the description of the observed recoil electron energy spectrum we 
have also considered the effect of departing from the SSM of Bahcall and Pinsonneault 1998 
by allowing arbitrary ^B and hep fluxes. Our results show that this additional freedom does 
not lead to a significant modification of the oscillation parameters. The best fit is obtained 
for ^B/^BssM = 0.61 and hep/hepssM = 12 for the SMA solution both for conversions 
into active or sterile neutrinos and ^B/^Bssm = 1-37 and hep/hepssM = 38 for the LMA 
solution. 
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TABLES 



Experiment 


Rate 


Ref. 


Units 


^BP98 


Homestake 
GALLEX ± SAGE 
Super-Kamiokande 


2.56 ±0.23 
72.3 ±5.6 
2.45 ± 0.08 


i 
ii 

i 


SNU 
SNU 
10^ cm-2 s"^ 


7.8 ±1.1 
130 ±7 
5.2 ±0.9 



TABLE L Measured rates for the Chlorine, Galhum and Super-Kamiokande experiments. 



Angular Range 


Dataj ± (Tj 


Day < cos 6* < 1 
Nl -0.2 <cos6l <0 
N2 -0.4 < cos 6* < -0.2 
N3 -O.6<cos0< -0.4 
N4 -O.8<cos0< -0.6 
N5 -1 <cose< -0.8 


0.463 ±0.0115 
0.512 ±0.026 
0.471 ± 0.025 
0.506 ±0.021 
0.484 ± 0.023 
0.478 ± 0.023 



TABLE IL Super-Kamiokande Cohaboration zenith angle distribution of events 
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Energy bin 


Dataj ± di^stat 


(^i,exp (%) 


Cri,cal (%) 


^i,uncorr (,/oj 


5.5 MeV <Ee<6 MeV 


0.472 ± 0.037 


1.3 


0.3 


4.0 


6 MeV <Ee< 6.5 MeV 


0.444 ± 0.025 


1.3 


0.3 


2.5 


6.5 MeV <Ee<7 MeV 


0.427 ±0.022 


1.3 


0.3 


1.7 


7 MeV <Ee< 7.5 MeV 


0.469 ± 0.022 


1.3 


0.5 


1.7 


7.5 MeV < ^e < 8 MeV 


0.516 ±0.022 


1.5 


0.7 


1.7 


8 MeV <Ee< 8.5 MeV 


0.488 ± 0.025 


1.8 


0.9 


1.7 


8.5 MeV <Ee<9 MeV 


0.444 ± 0.025 


2.2 


1.1 


1.7 


9 MeV <Ee< 9.5 MeV 


0.454 ±0.025 


2.5 


1.4 


1.7 


9.5 MeV < ^e < 10 MeV 


0.516 ±0.029 


2.9 


1.7 


1.7 


10 MeV <Ee< 10.5 MeV 


0.437 ±0.030 


3.3 


2.0 


1.7 


10.5 MeV < ^e < 11 MeV 


0.439 ± 0.032 


3.8 


2.3 


1.7 


11 MeV <Ee< 11.5 MeV 


0.476 ± 0.035 


4.3 


2.6 


1.7 


11.5 MeV <Ee<12 MeV 


0.481 ± 0.039 


4.8 


3.0 


1.7 


12. MeV <Ee< 12.5 MeV 


0.499 ± 0.044 


5.3 


3.4 


1.7 


12.5 MeV < ^e < 13 MeV 


0.538 ± 0.054 


6.0 


3.8 


1.7 


13 MeV <Ee< 13.5 MeV 


0.530 ± 0.069 


6.6 


4.3 


1.7 


13.5 MeV <Ee<U MeV 


0.689 ± 0.092 


7.3 


4.7 


1.7 


14 MeV <Ee<20 MeV 


0.612 ±0.077 


9.2 


5.8 


1.7 



TABLE III. Recoil energy spectrum of solar neutrinos from the 825-day Super-Kamiokande 
Collaboration data sample Q. Here (Ji^stat is the statistical error, ai^exp is the error due to correlated 
experimental errors, Ci^cai is the error due to the calculation of the expected spectrum, and (Tj, 
is due to uncorrelated systematic errors. 



'i,uncorr 



Period (month) 


Dataj ±ai (E > 11.5) 


0.0<t< 1.5 


0.588 ±0.057 


1.5 <t< 3.0 


0.588 ± 0.057 


3.0 < t < 4.5 


0.532 ± 0.069 


4.5 < t < 6.0 


0.392 ± 0.059 


6.0 < t < 7.5 


0.473 ± 0.059 


7.5 < t < 9.0 


0.521 ±0.065 


9.0 < t < 10.5 


0.548 ± 0.065 


10.5 < t < 12.0 


0.522 ± 0.058 



TABLE IV. Seasonal distribution of events given by the Super-Kamiokande Collaboration 
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FIG. 1. Theoretical recoil electron energy spectrum obtained by our calculation in Eqs. ^ 
([iq ) (solid histogram) compared with the Super-Kamiokande MC predictions (dotted histogram) 
Also shown are the data points from Super-Kamiokande data |P]. 
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FIG. 2. Allowed regions in Arn^ and sin from the measurements of the total event rates 
at Chlorine, Gallium and Super-Kamiokande (825-day data sample) for active-active transitions 
(a) and active-sterile (b). The darker (lighter) areas indicate the 99% (90%)CL regions. Global 
best-fit point is indicated by a star. Local best-fit points are indicated by a dot. 
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FIG. 3. Same as previous figure but allowing a free ^B flux normalization. 
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FIG. 4. Allowed regions in Arri^ and sin 6 from the measurements of the total event rates at 
Chlorine, Gallium and Super-Kamiokande (825-day data sample) combined with the zenith angle 
distribution observed in Super-Kamiokande for active-active (a) and active-sterile transitions (b). 
The darker (lighter) areas indicate 99% (90%)CL regions. Global best-fit point is indicated by 
a star. Local best-fit points are indicated by a dot. The shadowed area represents the region 
excluded by the zenith angle distribution data at 99% CL. 
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FIG. 5. Allowed regions in Am? and sin 9 from the measurements of the total event rates 
at the Chlorine, Gallium and Super-Kamiokande (825-day data sample) combined with the re- 
coil electron spectrum data observed in Super-Kamiokande for active-active oscillations (a) and 
active-sterile oscillations (b) . The darker (lighter) area indicate 99% (90%)CL regions. Global 
best-fit point is indicated by a star. Local best-fit points are indicated by a dot. The shadowed 
area represents the region excluded by the spectrum data at (99%)CL. 
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FIG. 6. Expected normalized recoil electron energy spectrum compared with the experimental 
data Q. The solid line represents the prediction for the best-fit SMA solution with free ^B and 
hep normalizations (/3 = 0.61, 7 = 12), while the dashed line gives the corresponding prediction 
for the best-fit LMA solution (/? = 1.37, 7 = 38). Finally, the dotted line represents the prediction 
for the best non-oscillation scheme with free ^B and /lep normalizations (/? = 0.45, 7 = 13.5). 
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FIG. 7. 95 and 99% CL regions obtained from the analysis of the seasonal dependence of the 
event rates observed in Super-Kamiokande from data with E^ > 11.5 for active-active oscillations 
(a) and active-sterile oscillations (b). The larger area represents the allowed region at 95% CL. 
The darker area is the excluded region at 99% CL. 
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FIG. 8. Allowed regions in Am? and sin 6 from the measurements of the total event rates at 
the Chlorine, Gallium and Super-Kamiokande (825-day data sample) combined with the zenith 
angle distribution observed in Super-Kamiokande, the recoil energy spectrum and the seasonal 
dependence of the event rates, for active-active oscillations (a) and active-sterile oscillations (b) . 
The darker (lighter) areas indicate 99% (90%)CL regions. Global best-fit point is indicated by a 
star. Local best-fit points are indicated by a dot. 
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FIG. 9. Same as previous figure allowing free ^B and hep neutrino flux normalizations. 
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